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variation which caused it. Two points appear at once on examination of the 
theoretical curves of figs. 2 and 3 : — 

{a) The galvanometer deflection does not give a perfectly faithful, but a 
slightly distorted, reproduction of the E.M.F. curve. 

(6) The galvanometer deflection lags behind the E.M.F. by an amount here 
equal to about one-fourth of the period. 

Both these deviations from a true representation of the E.M.F. variation 
are due to the effect of the inertia of the moving system of the galvanometer. 
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Section 1. — Note on the Direct Nmnerical Oalculation of Elliptic Integrals. 

The importance of realising rapid and accurate methods of calculating the 
elliptic integrals, now denoted by 

F (<^, h) = {^d^/A (^, k), E (<^, ^) = f A (<jb, k) d<j,, (1) 

Jo Jo 

where A (^, k) = y/' (1 — F sin^ ^), 

was first remarked by Euler* (1766), although it was not until several years 
later that Landenf (1775) discovered in geometrical form the transformation 
which forms the basis of existing methods of numerical calculation of the 
elliptic integrals. 

A method of successive transformations for the ultimate reduction of the 
algebraic forms of the elliptic integrals to elementary forms was published 
by Lagrange:!: in 1784-5. This memoir contains an exposition of the scales 
of arithmetico-geometrical means and discusses their use in calculating the 

*^ Enler, 'Novi Coinm. Acad. Sc. Petrop.,' vol. 10 (1766). 

t Landen, 'Phil. Trans. Roy. Soc.,' vol. 65, p. 283 (1775); also, 'Mathematical 
Memoirs,' London, 1780. 

J Lagrange, " Sur une nouvelle m6thode de calcul integral pour les differentielles 
afFectees d'un radical carr6 sous lequel la variable ne passe pas le quatrieme degr6," 
'Mem. de FAcad. Koy. des Sc. de Turin,' vol. 2 (1784-5) : ' CEuvres,' vol. 2, pp. 253-312 
(Gauthier-Villars, Paris, 1868). 
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elliptic integrals of the first and second kinds in a manner }3ractically 
identical with that formulated about the same time by Legendre.*' 

This method, discovered independently by Gauss and called by him the 
"algorithm of the arithmetico-geometrical mean," originated in connection 
with the evaluation of complete elliptic integrals arising from a problem in 
attractions occurring in the planetary theory.f In a note he mentions that 
the results were obtained by him as part of a more comprehensive theory 
independently of the results of Lagrange and Legendre, with which they are 
closely connected. 

Much use was made of the scales of arithmetico-geometrical means in the 
numerical calculation of elliptical integrals by Jacobi,| who, however, adapted 
his formulse to computation by logarithms. The use of modern calculating 
machines has greatly extended the scope of direct numerical calculation of 
elliptic integrals by the methods just referred to. In many practical applica- 
tions it is advantageous to have available a direct method of computation 
independent of auxiliary tables, the use of which often requires tedious and 
somewhat uncertain interpolations when seven or more significant figures are 
desired. Many formulae involving elliptic functions arising in mathematical 
physics, when expressed directly in terms of quantities arising from the 
arithmetico-geometrical scales, give rise to extremely convenient expressions 
for direct numerical computation when a calculating machine is available. 

Section 2. — The Scale of Arithmetico-geometrical Means. 

In forming the scale of arithmetico-geometrical means, we start with the 
positive numbers, % = 1, 5o = k\ cq z=: k, k being the modulus of formula (1), 
and k' the complementary modulus, k'^ = 1— F. We calculate successively 



^0 = 1, 


ho- 


^k\ 


Cq = 


^k 


^1 == i(%+&o) 


h-- 


= v^C^o^o), 


Ci = 


= |(%-^o) 


a2 = |- (^1 + h) 


h-- 


= \/(«i^iX 


^2 = 


= i(%— ^i) 



N 



y* (ijj 



% = i (^»-l + &»-lX hn = \/ iftn-l &»-l), C^ = I (^n-l — &n-l) 



••••• •>«•*• 



* Legendre, " M^moire sur les integrations par arcs d*ellipse," — " Second memoire, 
etc.," * M6ni. de FAcad. des Sciences de Paris, annie 1786,' pp. 616-643 (Paris, 1788) ; 
also pp. 644-683 ; *Trait6 des Fonctions Elliptiques,' vol. 1, pp. 79 et seq. (Paris, 1825). 

t Gauss, " Determinatio Attractionis, etc.," ^Comm. Gott. Soc. Eeg. Scient.,' vol. 4, 
1818 ; * Werke,' vol. 3, pp. 331-355 (Gottingen, 1866). The note mentioned above is 
dated February 9, 1818 j ' Werke,* vol. 3, pp. 357-360. 

J Jacobi, * Fundamenta Nova,' §§ 38 and 52 (1829) ; * Ges. Werke,' vol. 1, pp. 154 
and 203 (Berlin, 1881) ; also, " Numerische Berechnung der EUipfcischen Functionen," 
* Crelle,' vol. 26, pp. 93-114 ; * Ges. Werke,' vol. 1, pp. 345-368 (1881). 
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The following relations are easily verified — 

The array of numbers (2) will be briefly referred to as the A.G.M. scale 

(aQ = 1, 60 = k'). The successive a's and 5's converge with extraordinary 

rapidity to the same limit, which Gauss denoted by M (aa, ho). It is easily 

seen that 

M.{€ao, eho) = €M{ao, ho), (4) 

The limit M.(aQ, ho) will be denoted by a^j as long as by so doing no 
ambiguity is involved. 

The following expressions for the complete elliptic integrals in terms of the 
A.G.M, scale (ao = 1, &o = ^') are well known. 

K = I tt/^w, (5) 

(K-E)/K= |(co^ + 2ciH4c2H... + 2X^+...). (6) 

It does not appear to have been remarked hitherto that 2""" log (^n/c„) also 
converges rapidly to a definite limit which enables us to compute the comple- 
mentary integrals K', E', and the nome q given by log^ = — ttK/K'. 

It may be proved* that 

Iwan/a^' == iw W/K- 2-'' log (a^/cn). (7) 

More convenient for calculation is the following formula, derived from (7) 
and the relations (3) 



00 



log (I/2) = TT K'/K = log (4ai/ci)-2 (|)»-i log («„/a„+i). (8) 

Making use of Legendre's relation, 

E/K + E'/K'-l = i7r/(KK'), (9) 

we easily derive the series 

K'-E' = Ja„(log^^-log^-i log^^-... 

(1^1 C^2_^^3_2C22-~4C32-. ..)-«„. (10) 

"^ The proof will be given in a memoir on the numerical calculation of elliptic functions 
which the writer has in preparation. The writer has been unable to find this formula 
explicitly stated by Gauss or by any subsequent writer on elliptic functions. It 
furnishes the key, however, to many remarkable formuUe found among Gauss's 
unpublished papers, stated without proofs (* Werke,' vol. 3 (1866)); Nachlass, 
' Arithmetisch-Geometrisches Mittel,' pp. 361-402 ; also, * Zur Theorie der Neuen 
Transscendenten,' pp. 446 et seq.). The series (8) is, however, given without proof by 
Gauss {loc. cit, p. 377), who may have been in possession of (7), but for some reason 
disposed to withhold the publication of this key-formula. 
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Section 3.- — Amplication to Mutual Inductance Formmlcefor Coaxial Circles, 

Let A be the radius of the larger circle^ a that of the smaller, and d the 
distance between their planes. We also denote by ri and rg the greatest and 
least distances respectively of a point on one circle to the circumference of 
the other. Maxwell* derived two well-known formulae for M, the coefficient 
of mutual induction of the two circles ; 

M = 87r(K//^)v/(A^)[l--E/K~|Fj "^ 

where k^ = 4:Aa/ {{a + Af + d^}, // = ^2/^1 K (11) 

and n^ = (A-i-ay + d^, ri = {^A—Gbf-^-d^ ^ 

Landen's transformation applied to (11) gives rise to Maxwell's second 

formula 

M = Stt v/(A^/^i)(Ki-~Ei), (12) 

where Ki and Ei are the complete elliptic integrals to modulus 

fti = (n — T2)l{r^ + r2), (13) 

If we compute the A.G.M. scale {a^ = n, 60 = ^'2), it is readily proved from 
(4), (5), and (6) that both (11) and (12) give rise to 

M = {27rVa4[ci2 + 2c2H4c3H... + 2^-ic„H...], (14) 

a most convenient and elegant formula for the direct computation of M. 

This series is always extremely convergent, since from (7) the successive 
values of c tend to converge like the series 

X x^ x^ x^ x^^...x^''.,. where a:; <1. 

When the circles a>ve very close together, ¥ is small and it is more convenient 
to use the complementary A.G.M. scale, which is then extremely convergent. 
The first few terms are as follows : — 
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an 


bn' 


^n 





n+^2 


Ti^T2 


2v^(nr2) 
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2^{Aa) 


^3 


2 


|n+v^(Aa) 


(4Aari2)l 


In— v^(Aa) 



\ 



V 



J 



(15) 



Making use of (10) and (12), we easily obtain the following formula for M 
in terms of the complementary A.G.M. scale (denoted by accented letters) : — 



M = 27raJ 



r (,,._2c/--4c3-- ^) (^i^g 4^_ ^ iog£/_ jiog^ 



an \ rg ^3 ^4 

(16) 
Logarithms to the base 10 may be used in computing the second series, 

* Maxwell, * Electricity and Magnetism,' 3rd ed., see. 701 (1904). 
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which may finally be multiplied by the factor fi == logglO = 2'30258509 to 
convert to logarithms to base e, 

For accuracy the first term of the above formula should exceed twice the 
second term. This will be the case if Tijr^ > 13, about. Hence formula (16) 
is suitable for computation when the circles are so close together that 
^'2/^1 <C O'l- It may, without serious loss of accuracy, be extended to cases 
where r2/ri<:^0'2y in which case seven-figure accuracy only requires the 
complementary A.G.M. scale (15) to be computed as far as n = 2. 

The formula (14) is, however, so much simpler than (16) that it may be 
used even when the circles are very close together. It has the advantage of 
depending only on figures derived by direct computation on a calculating 
machine, and does not depend on the use of logarithmic tables. 

By dividing up coils of finite cross-section into a number of coils of 
sufficiently small cross-section, we may make use of a formula of quadratures 
for computing the mutual induction, depending on the calculation of the 
coefficients for the elementary coils taken in pairs. For this purpose the 
simple formula (14) is specially applicable, while (16) may be useful in 
some cases. Extremely convergent formulae for the numerical computation 
of the induction coefficients of coaxial single-layer coils may be obtained by 
the use of A.G.M. scales and trigonometrical recurrence-formulae. These 
depend on a number of new formulae for the computation of the incomplete 
elliptic integrals of the first, second, and third kinds, which it is hoped to 
deal with in a future paper. 

Section 4. — Numerical Ilhtstrations. 

To illustrate the rapid convergence of the A.G.M. scales, Gauss* works out 
an example, starting with % = 1, &o == 0*2, and finds that % and ^5 are 
represented by the same number to fifteen significant figures. 

The usual arithmetical process of extracting square roots is easily adapted 
to calculating machines of modern construction.f In carrying out the 
computations, the work is considerably shortened by keeping in mind the 
rule that if the first p out of the n digits of the number required in the 
square root have been obtained by the usual process, the next (|?-"1) digits 
can be obtained by division only, with a possible error of 1 in the last digit.^ 

The following numerical examples, taken from the collection of formulae 

* Gauss, ' Werke,' vol. 3, p. 363 (1866). 

t A description of modern calculating machines is given by Horsburgh, E. M., 
'Modern Instruments and Methods of Calculation' (London, G. Bell, and Sons, 1914) ; 
also by d'Ocagne, * Le Calcul Simplifi6' (Gautliier*-Villars, Paris, 1905). 

X Chrystal's 'Algebra,' Part I, 5th ed., p. 210 (1904) 
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published by the Washington Bureau of Standards,* illustrate the advantages 
of formulae (14) and (16) in numerical computation. 

\ Example 1. ' Bulletin 169/ p. 32 (1916). Equal circles at a distance, 
A = a = 25 cm., <i = 50 cm., ri = 50^^ 2 cm., r^ = 50 cm. 



n. 


an. 


in. 


Gn. 


C'n. 



1 
2 
3 


70 -71067812 
60 -35533909 
59 -90784740 
59-90701174 


50 -00000000 
59 -46035575 
59-90617607 
59-90701173 


10 -35533906 
-44749167 
-00083567 


107 -23304776i 
-200248795 
-000000698 



To an order of accuracy of ten significant figures, we need take only three 
terms of (14), which then gives, with 27r^ = 197392088022, 



TT 



M = ^^(ci2 + 2c22 + 4c32) 



_ 27r2x 107-633548143 



= 35-464981878 cm. 



a^ • ■ 59-9U7011736 

It is readily seen that little additional labour is involved in carrying out 
the computations to ten significant figures. Values quoted in the above- 
mentioned * Bulletin 169,' computed by three different formulae expressed in 
^-series, derived by Nagaoka, give for M the three values 

35-464975 cm., 35-464981 cm., 35*46481 cm. 

In carrying out the computation, there is no need to write down more 
figures than are set down above, all the processes being carried out on the 
calculating machine. 

Example 2. 'Bulletin 169,' pp. 21, 22, 28 (1916). Equal circles close together, 
A = a = 25 cm., cZ = 1 cm., ^ n = 50-010000. r^ = I'OOOOOOOO. 
From (15) we write down the complementary A.G.M. scale — 



n. 


dn . 


hn\ 


G W 


1 
2 


50 -010000 
50 -OOoOOO 


50 -000000 
50 -005000 


1 -000000 
-005000 



Formula (16) then gives to eight significant figures. 



M = 27r 



\ ^3 / ^2 



TT 



2501-0000 
L50-005000 



X log (200-02000) -100-01000 



or 



M 



TT 



1000-4000 

Lio-001000 



X 5-2984174-- 200-0200 



1036-66486 cm. 



"^ Bulletin 169, * Formulas and Tables for the Calculation of Mutual and Self- 
Inductance,' by Kosa and Grover, 3rd ed., 1916. 

VOL. C. — A. F 
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Values derived from Legendre's tables and two series formulae are given as 

1036-6664 cm, 1036-6652 cm, 1036-663 cm. 

A comparison of the numerical work involved emphasises in marked 
manner the advantage of the new formula (16). 

It is interesting to check this result by formula (14), although the 
present example is not favourable to its use. 

The computation proceeds as follows : — 



n. 


(ln» 


hn. 


Cn. 


Cn'^. 



1 
2 
3 
4 
5 


50 -OlOOOo 
25 -50500o 
16 -288387 
14 -859193 
14-82475 
14 -82473 


1 -OOOOOo 
7 -071774 

13 -43000o 

14 -79030 
14 -82471 
14 -82473 


24 -50500 
9 -21661 
1 -42919 
-03445 
-00002 


600 -4950 

84 -9459 

2 -0426 

-0012 

-0000 



Even ill this unfavourable case, four terms of (14) are more than sufficient 
to give an accuracy to seven significant figures, and 



M 



2t' / 2 , o 2 , ^ 2 , e 2N 27r2x 778-5666 



= 1036*666 cm. 



The labour of computation by the above formula is still less than that 
involved in the use of Legendre's Tables or of any of the series formulae. 

Section Z,— Summary and Conclusions, 

The use of the modern calculating machine for the extraction of square 
roots makes the scales of arithmetico-geometriccd means extremely convenient 
for the direct and rapid computation of the complete elliptic integrals. 

This procedure is applied to Maxwell's formulae for the coefficient of 
mutual induction of coaxial circles, with the result that two new and 
elegant formulae are obtained for numerical computation. 

Numerical examples, taken from the well-known collection of formulae 
and tables for the calculation of self and mutual induction, issued by the 
Bureau of Standards, illustrate the advantage of the new formulae in 
economy of labour. 



